The role of entanglement in work extraction 
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We consider reversible work extraction from identical quantum batteries. From an ensemble of 
individually passive states, work can be produced only via global unitary (and thus entangling) 
operations. However, we show here that there always exists a method to extract all possible work 
without creating any entanglement, at the price of generically requiring more operations (i.e. ad- 
ditional time). We then study faster methods to extract work and provide a quantitative relation 
between the amount of generated multipartite entanglement and extractable work. Our results 
suggest a general relation between entanglement generation and the power of work extraction. 
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Energy storage and its subsequent extraction has al- 
ways been a central topic of thermodynamics, due to 
its obvious fundamental and practical importance. En- 
tanglement is a key feature of quantum mechanical sys- 
tems. It is one of the cornerstones of quantum infor- 
mation theory and it proved to be important in differ- 
ent physical phenomena as, to name a couple, quantum 
phase transitions [U [2] or fractional quantum Hall effect 
[5]. In quantum thermodynamics, entanglement is also 
connected to work extraction from multipartite systems. 
Indeed, global (and thus entangling) unitary operations 
are capable of extracting more work than local opera- 
tions from a set of quantum batteries (see, for instance, 
[3]). However, the fact that an entangling operation is 
needed for a process, does not imply that entanglement 
is generated during its execution. The scope of this work 
is to clarify these connections. We first show that op- 
timal work extraction can be achieved without any en- 
tanglement generation. The corresponding dynamics is 
slow, in the sense that it requires many different oper- 
ations. We then consider faster methods and establish 
a link between entanglement generation and extractable 
work. Our results point to a connection between entan- 
glement generation and work power. 

A classical result on work extraction is the no-go the- 
orem called " Thomson's formulation of the second law 
of thermodynamics: no work can be extracted from a 
thermal state via cyclic hamiltonian processes [5]. It is 
extended to finite quantum systems by introducing the 
concept of passive states, those not capable of giving 
out energy during any cyclic hamiltonian process [6l [7] . 
Thermal states are passive, while not all passive states 
are thermal. Indeed several copies of states of the lat- 
ter type may contain population inversions, allowing for 
some work extraction. On the other hand, thermal states 
do not share this ability of " activation" since any com- 
bination (tensor product) of thermal states is a thermal, 
thus passive, state. The converse is also true: if any com- 
bination of a passive state is passive then it is a thermal 



state. One can, however, get some work out of a set of 
locally thermal states by initially correlating them, the 
microcanonical state being a prominent example [8111 lj . 

We study the problem of work extraction from N non- 
interacting finite-level batteries. For simplicity and with 
no loss of generality we assume all batteries to be identi- 
cal d- level systems with hamiltonian H = J\ =1 £k\k){k\, 
£i < ■•• < £d- The hamiltonian of the ensemble then 
reads as h = H (g) • • • <g) 1 + ... + 1 (g) • • • <g) H. We take 
the eigenbasis of ho as the standard basis. 

Typically one wishes to store (even if not for long time) 
the energy before extracting it. For that, one prepares 
the ensemble in an active state which is diagonal in the 
total hamiltonian so that it remains unchanged: 

n = diag(pi, p d N) (1) 

d N 

with all p a > and ^2 a=1 p a = 1- Note that the initial 
state is separable, as it is diagonal in the product energy 
eigenbasis, but may display classical correlations. Rele- 
vant examples of states having this form are the micro- 
canonical state [8HTT] or a set of passive-but-not-thermal 
states HH[7]. 

Now, to extract work, at the moment t = exter- 
nal control fields with a time dependent potential V(t) 
are turned on. In the presence of these fields the sys- 
tem — made up of N batteries — undergoes a unitary evo- 
lution driven by the hamiltonian 

h(t) = h + V(t). (2) 

Then, at the moment t = r, external fields are turned 
off: V(t) — 0; providing, thus, the cyclicity of the work 
extraction process and returning the system its original 
structure. The overall unitary evolution operator of the 
system is given by the time-ordered exponent U(t) — 
exrj (— i fl dt [h + V(t)}) , rendering the final state to be: 

fl( T ) = U{T)nu\r). (3) 
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The work extracted during the process is given by the 
difference of the initial and final energies of the system: 

W^ = tr(n/i )-tr(fi(r)/io). (4) 

As proven in 8J, the unitary operator minimizing 
tr(f2(r)/io) (hence maximizing W) is the one that re- 
sults in a permutation of the elements of VL so that the 
largest element of f2 is matched with the smallest one 
of ho, the one-but-largest with the one-but-smallest, etc. 
This shows that the state is passive if its elements are 
reversely ordered with respect to the hamiltonian ho- 

Now note that after maximal work extraction the final 
state of the ensemble is separable because is diagonal in 
the eigenbasis of the non-interacting hamiltonian ho, as 
the initial state. On the other hand, permutation op- 
erations are not local and can have maximal entangling 
power [12] . so the following natural questions arise: does 
the state of the system get entangled during the process? 
If yes, how entangled does it become? Is there any way 
to bypass the entanglement creation, so that the state 
remains classically correlated all the time? 

We start by answering the last question and provide a 
protocol that attains maximal work extraction with no 
entanglement. Then, we consider faster protocols and 
provide lower bounds to the entanglement they gener- 
ate. These bounds are applied to two relevant scenarios, 
namely maximal work extraction from the microcanoni- 
cal bath and from copies of passive states. 

Bypassing entanglement. Indirect paths. — To extract 
maximal work from a state of the form ([!]) wc need to 
reorder its entries accordingly. This reordering can be 
done in elementary steps of transpositions. E.g., may 
be such that the population of the lowest energy level 
(pi) is not its maximal element, namely p^i- Then 
one needs to transpose p\ with p^. If needed, one will 
have to transpose the population of the one-but-lowest 
energy level with the one-but-largest element of f2, etc. 
After some number of such steps the state will be ordered 
properly, yielding maximal work. 

The transposition of the population of any energy level 
by the population of some other level can be done without 
creating entanglement in the meantime. Indeed, suppose 
we need to transpose p a by pp. The eigenstates of ho 
that correspond to them are \a) = Iz"^...*^) and \(3) = 
| if %2 ...i^ N ), respectively. We then divide this action in 
2N — 1 transposition steps. First 

Kl*2— i%) Kl*2"-*jv) ^ 1*1*2 1*1*2— *'jV") (5) 

and the -/V— 1 steps back from \if i^ ...i^_ 1 i%) to \a). On 
each step one exchanges the populations between states 
involved in it. Here we stress that all steps in ^ involve 
only the corresponding states, while the populations of 
the other states are kept unchanged. Thus, one cannot 
perform any if these steps with a local unitary operation. 



In general, the unitary transposing the populations of 
two basis states, say \a) and |/3) reads as 

U aP = £ | 7 )<7l + H</3| + |/3)R (6) 

If the control potential V(t) ^ generating U al3 couples 
only to | a) and \/3), the evolution operator at some in- 
termediate moment t of the process is 

U aP {t)= £ l7><7l+u°*(*), (?) 

where u a ^(t) lives in the linear span of \a) and |/3) and is 
unitary. It depends on t and the concrete form of V(t). 

Now for, e.g., the first step in ^ we need to perform 
the transposition unitary U aa between \a) — \i"i% 
and \a!) = |if i% ...ijv). According to |7| the state of the 
ensemble at an intermediate moment t is 

Q(t) = U aa ' ' (t)VLU aa '\t) 

={p a +Po l ')pi{t)®\i^...^ f )^...i%\ + Y J P7l7><7l-(8) 

Quite straightforwardly, pi(t) > and tr[pi(i)] = 1, so 
pH) means that the state of the ensemble is separable 
during the whole process of population exchange between 
\a) and \a'). Notice that although U aa (t) is global and 
thus has entangling power, there are states which do not 
become entangled under its action. 

By the same reasoning, one may stay separable also 
during the rest of transpositions in chain ([5]). So any re- 
placement in the global state of the battery can be made 
without creating entanglement between its constituents, 
which proves that one can extract maximal work from 
the battery and stay separable during the whole process. 

Finally we note that the results of this section also hold 
for a generic product initial state ®iPi, where each pi may 
have quantum coherences, that is, off-diagonal elements 
in the energy eigenbasis. One first extracts maximal work 
from individual systems via local, thus non-entangling, 
unitaries, and then applies the protocol described above 
to the resulting diagonal state. 

The previous non-entangling protocol requires 2N — 1 
global operations in order to perform the desired ex- 
change of the populations \a) and |/3). However, this 
exchange can be performed in one step by the unitary 
We term such evolutions by direct paths. Now, the 
natural question is whether entanglement is generated by 
these direct paths, which allow one to extract work faster 
and thus get more power. 

Direct paths. — Consider the population exchange of 
| a) and A relevant example of a direct path 

would be the time independent hamiltonian H = 
7p (|a)(/3| + |/3) (a |) which generates the desired inter- 
change at t = t. This corresponds to a constant field 
acting only on \a), |/3) for t € (0,r). More generally, we 
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will consider the evolution of Q(t) — U a @Q,U a ^ where 
U a0 is found from Q. 

We will take all indices of \a) and |/3) to be different. 
Suppose \a) and \j3) would only have n different indices. 
By proceeding as in we notice that time evolution 
only creates coherences in the corresponding n batteries, 
so only they can get entangled. Therefore all results we 
will derive hold in the general case: one just needs to 
replace N by n. 

In order to measure entanglement of tt(t) in a direct 
path we use a recently proposed computable measure of 
genuine multipartite entanglement for mixed states [T3] 
which employs the notion of /-separability. A pure state 
is /-separable if it can be written as a tensor product 
of / terms. Then, a mixed state is called /-separable if 
it can be decomposed into pure states that are at least 
/-separable [22] . 

Now, to introduce the mentioned measure, consider a 
pure TV-partite state ip. We define the entropy vector 5 
53 as 

S fc (V0 = ^2(l-tr(p2)), k=l,...,2 N ~ l -l, (9) 



where the index fc runs over all possible bipartitions 
r = {{"fk\fk)} of {1, 2, ...,N}; and p k is the corresponding 
reduced state: pk — Tr 7 - fe (|V')(V'l)- If we se t tne entries 
S k in non-increasing order, then it holds that (i) if the 
last 2 i_1 — 1 entries are zero, then the state is /-separable, 
(ii) the last entry is a measure of genuine iV-partite en- 
tanglement. 

This measure can be extended to mixed states p via a 
convex roof construction |15j : 

E k (p) = inf VftSfctyi), fc = l,...,2 Ar - 1 -l, (10) 

{pi-,\-4>)}^ 

where the infimum is taken over all possible decomposi- 
tions {pi\ \tjj)} of p, and Sk are the entries of the entropy 
vector arranged in non-increasing order. Then, the state- 
ments (i) and (ii) above also hold for E k (E k plays the 
role of Sk in mixed states) |15j . 

From Q we notice that the state Q(t) has only two 
non-zero off-diagonal elements and four time dependent 
elements. For such a simple form the lower bounds from 
Ref. [15] for Ek read as 



E k (Q(t)) > 2 



(a\n\a) 
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<fc) 



\a)^\P lt nPi\fi) 



(11) 



where {Tj '} is the set of all subsets of V containing 
k elements, Y = {(7i|7i)} is the set of all possible bi- 



partitions, and Py = 



Dm 



The right hand 



side of (11) can be made time-independent by using 
max t {\uf^{t)\\u"^ 2 (t)\} = 1/2 — a consequence of its uni- 
tarity. 

Therefore, we find conditions on /-separability for the 
exchange of the population of two states that depend 
only on the initial state ft. Indeed, if -E 7 (z) > 0, with 
7(Z) = 2 Ar ~ 1 — 2' -1 + 1, then we know the state is at 
most /-separable during any direct path. Furthermore, 
the bounds are tight for the first and the last entry [IB] , 
and thus E\ and E 2 n~i_ 1 can be used to detect entan- 
glement and measure genuine N— partite entanglement 
respectively because of statements (i) and (ii). 

Microcanonical bath — As a first application of our 
derivations we consider work extraction from the micro- 
canonical state of width A around some energy E of the 
total system 
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E 

ni—l 



(12) 



with \Eq ) being states with energy in the interval [Eq- 



A/2, E + A/2] [33]. This state is active [gHTT] but is not 
capable of giving macroscopic work in thermodynamic 
(N > 1) limit [IT] . 

In this setting, work is extracted by performing 
exchanges of the populations (all equal to 0) of the Na 
states |eo"^) (m = 1,...,Na) with the lowest energy by 
the populations (all equal to 1//Va) of the states |i?o"^)- 
These exchanges do not overlap and can be done in Na 
successive steps. So the bounds ( [TT| ) can be used straight- 
forwardly for each m. If as a first exchange we take the 
one between states with energy and Eo — A/2, then if 
there are ri\ different indices in those sates [23, formula 



(111 shows that the corresponding bath particles are gen- 



uinely n!-partite entangled during the whole transposi- 
tion. [35] 

A set of passive states — Another important scenario is 
the product of passive-but-not-thermal states [3 El E] : 



* L ap — u p 



(13) 



with 



diag(pi, ...jPd). The amount of work is 



bounded by the free energy difference [3: 

W = iVtr [H(a p - a th )] , (14) 




where a th is thermal and has the same entropy as a p . 

For such a product form our results become partic- 
ularly simple. Consider the direct path exchange of 
the populations of \a) and \/3) given, respectively, by 
p a = Y[i=iPi a an d P/3 — rii=i?'i' 3 - Here the expressions 



(111 are simplified to 



E k (fl(t)) > \pa~Pp\- ZkyfpaJ 



(15) 



where we also maximized the right hand side over t. 
Therefore, the state is at most Z-separable if: 



Pa 



> l + 2 7 + 2 V7T72, 7 = 2 W " 1 - 2'- 1 + 1; (16) 



for the Pa > P/3 case. 

It was recently shown [4] that in the thermodynamic 
limit (N 3> 1) the ultimate bound (14) for extractable 
work is asymptotically reachable via a global unitary (de- 
note it by Un). In the construction of Un the key idea 
is that there are the same number of typical states in 
<S) op and cr because their entropies are the same. 
Hence the probabilities of these typical states coincide 
also. Since the overall probability of nontypical states is 
exponentially small, upon transposing in (g) N o~ p its typi- 
cal states with the states that are typical in ® N a i , we 
arrive at an energy approximately equal to the one of 
® N a th . These e NS ^ transpositions can be made suc- 
cessively and without overlaps. 

Typical states of ® N a p are such \ii, ...ijv)s that contain 
(poN+N ) zeroes, (piN+Ni) ones, etc., where Yl^i = 
and all |JVj| < sN with arbitrarily small e > 0. Their 
populations are 



IK Pfc IK 



(17) 



k=l 



k=l 



N a th consist of 



Analogously, the typical states of 
(pjjfW + N' k ) number of \k)s (k = 1, d), so their popu- 
lations in (S) N a v read as 



n 

k=l 



Pk 



(18) 



k = l 



with J2 N i=° and l-^il < eN - 

If we choose to transpose states with {jVtlLj = 
{N' k }f =1 , then, for all transpositions, formula (16) will 
imply the same condition 



In 



1 + 2 7 + 2^7 + Y 



N 



(19) 



that o~ p must satisfy to be at most /-separable during the 
proc ess. Here S( \ \ ) is the conditional entropy; 7 is from 
(16) and S(a p ) — S(<r th ) is employed. 



P, 
0.2 




FIG. 1: Contour plot of the work W released by four three- 
level systems initially in the state ® A o p on the exchange of 
populations pf and P0P2 of levels 1 1 1 1 1} and 1 0222} . Transpo- 
sition is done following a direct path. The levels of each sys- 
tem are {0, e, e}. Lighter regions correspond to more work ex- 
traction. The dashed lines separate regions of fc-separability, 
the left side (SEP) being fully separable and the right most 
region (GME) genuinely multipartite entangled. Notice that 
there are regions where work is extracted without generat- 
ing any entanglement during the process. However, when 
W is larger than a certain threshold, entanglement is always 
present. Furthermore, as we increase W, genuine i-partite 
entanglement starts appearing with I progressively larger. 



Condition ( 19 1 has a simple interpretation — the 
greater the difference between a p and cr th , the more en- 
tanglement we need. In the N ^ 1 limit, the con- 
dition for entanglement being present is 5 f (cr th ||o'p) > 
ln[3 + 2\/2]/-/V so basically all states get entangled, while 
the condition for genuine iV-partite entanglement to ap- 
pear is 5 , (cr th ||CTp) > ln[2] which tells that this entangle- 
ment does not have to be iV-partite. Remarkably, the 
latter condition is independent on TV. 

On the other hand, the extracted work is governed by 



the difference of a p and a th (14). So the farther a th is 



from a p the more entanglement is generated and more 
work is extracted. This situation is illustrated in Fig. [l] 
for four three-level systems. 

Entanglement and power — These two exemplary cases 
show that entanglement is widely present during direct 
routes. Furthermore, we have shown that the amount of 
entanglement is directly connected to the amount of work 
for the case of identical batteries. This statement is not 
only applicable to this case. Indeed, consider the work 
extracted in a two-state exchange, which is simply pro- 
portional to the difference of their populations. Notice 
then the entanglement generated is directly related to 
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this difference because of the first term in ( |11[ ) , although 
in this case the exact value of entanglement entropy also 
depends on the other diagonal elements (because of the 
subtracted term). 

On the other hand, the amount of entanglement gen- 
crated may be reduced by following a combination of di- 
rect and indirect paths. Indeed, one simply performs 
N — I populations exchanges via indirect paths as in ^ , 
followed by a direct path which will produce at most 
I— partite entanglement. Alternatively, for the case of 
identical batteries, where the amount of multipartite en- 
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tanglement generated is given by condition (19), one can 



arbitrarily reduce the amount of entanglement generated 
by performing K extrasteps involving the direct exchange 
of the populations of states pk (k = 1, K) which satisfy 
S{p\\\(Jp) < S(p2\\a p ) < ... < S(a th \\(Tp). In both cases, 
if we assume that all global transposition operations are 
equally time consuming, reducing entanglement produc- 
tion comes at the expense of increasing the time of the 
process (recall that the same permutation takes 2A — 1 
times more operations to perform via indirect paths than 
by direct exchanges). 

Our analysis thus suggests that entanglement plays no 
role for the amount of work one extracts but rather it is 
crucial for the power — the larger the power output the 
more entanglement is created during the process. This 
strictly holds for the protocols considered in this work. 

Conclusions & Outlook. — We studied the role of entan- 
glement generation in work extraction from an ensemble 
of N non-interacting batteries. We introduce a protocol 
for maximal work extraction such that no entanglement 
is created in the ensemble during the runtime. We then 
consider direct paths — defined as evolutions where the 
controlled potential is only coupled to the states from 
which one wishes to extract energy — and show that a 
higher amount of work extraction requires a higher level 
of entanglement (see, for instance, Fig. [lj. For the con- 
sidered protocols, the power at which work is extracted 
is connected to the entanglement production. 

For future investigations, it will be interesting to study 
usable work extraction versus the entangling power of the 
required unitaries. 

Also, the task of permuting different diagonal elements 
in a density matrix is essential in other scenarios in quan- 
tum thermodynamics. Indeed, recently there have been 
proposals of resource theories for usable work extraction 
via energy preserving unitaries Ref. |191l21j . There, work 
is also studied as a swap of states from an initial diagonal 
state, so our results regarding multipartite entanglement 
generation straightforwardly apply. Furthermore, simi- 
lar considerations hold for other problems like dynamical 
cooling of spins in NMR physics or constructing optimal 
heat engines (see, e.g., [T71 [TBI HO]), etc-, so our results 
have implications also beyond the problem of work ex- 
traction. 

We thank D. Cavalcanti, A. E. Allahverdyan, J. B. 
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one can connect the notion of /-separability to multipar- 
tite entanglement by noting that there is at least (N/l)- 
partite entanglement and at most (N — I + l)-partite en- 
tanglement in an /-separable state 

Although A must be much smaller than E, it must be 
large enough for In (-/Va) be oc N. 

For the bath consisting of N spins with gaps S it is easy 
to see that n\ = E ~^/ 2 with pi being the population of 
the excited level of one bath particle. 



Indeed, notice that the second term in (111 vanishes be- 
cause, due to Eo 3> A, at least one of the diagonal ele- 
ments under the square root will lay outside the energy 
range [E - A/2, E Q + A/2]. 



